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Executive  Summary 


The  proposed  research  has  been  focused  on  the  modeling  of  quantum  processes 
in  intersubband  lasers  of  Si/ZnS  superlattices,  which  could  potentially  offer  the  great 
opportunity  for  monolithic  integration  of  Si-based  photonic  devices  with  advanced  Si 
electronics.  There  are  many  new  important  and  interesting  quantum  effects  in  this  new 
mixed  polar-nonpolar  material  system,  arising  from  the  fact  that  not  only  the  carriers 
are  confined,  but  also  the  polar  and  nonpolar  optical  phonons.  So  far,  there  has  been 
little  reported  investigation  on  this  material  system.  During  the  period  of  this  research 
contract,  we  have  calculated  the  valence  subband  engeries  and  wavefunctions  in  Si/ZnS 
quantum  well  structures.  We  have  also  established  models  for  the  confined  optical 
phonons  in  Si/ZnS  superlattices.  Optical  phonon  modes  in  Si  and  ZnS  layers  are  totally 
confined  within  their  respective  layers  since  both  layers  can  be  treated  as  infinitely  rigid 
with  respect  to  the  other  layer.  Since  there  are  no  associated  electric  fields  with  nonpolar 
optical  phonons  in  Si  layers,  only  a  mechanical  boundary  condition  needs  to  be  satisfied 
for  these  nonpolar  optical  modes  at  the  Si-ZnS  interface.  The  optical  phonons  in  Si 
layers  can  be  described  by  guided  modes  consisting  of  an  uncoupled  s-TO  mode  and  a 
hybrid  of  LO  and  p-TO  modes  with  no  interface  modes.  In  ZnS  layers,  a  continuum 
model  hybridizing  the  LO,  TO  and  IP  modes  is  necessary  to  satisfy  both  the  mechanical 
and  electrostatic  boundary  conditions  at  the  heterointerface.  A  numerical  procedure  is 
provided  to  determine  the  common  frequency  between  LO,  TO,  and  IP  modes.  This  is 
a  new  procedure  for  obtaining  the  eigen  modes  of  a  mixed  polar-nonpolax  heterosystem. 
Analytical  expressions  are  obtained  for  the  ionic  displacement  and  associated  electric 
field  as  well  as  scalar  and  vector  potentials.  The  established  optical-phonon  model  has 
then  been  used  to  calculate  the  phonon  scattering  process  of  heavy  holes  in  Si/ZnS 
heterostructures.  Our  results  indicate  that  contributions  to  the  intersubband  scattering 
rate  from  Si  or  ZnS  confined  optical  phonons  depend  strongly  on  the  distribution  of 
envelope  wave  functions  over  the  respective  layers  within  which  different  types  of  optical 
phonons  are  confined.  This  work  established  a  theoretical  foundation  for  studying  the 
lasing  lifetimes  of  Si/ZnS  intersubband  lasers. 
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I.  INTRODUCTION 


The  demonstration  of  the  InGaAs/AlInAs  intersubband  quantum  cascade  laser 
at  A  =  4.2//m[l,  2]  has  spurred  interest  in  the  use  of  silicon  as  the  lasing  material 
because  of  its  integrability  with  advanced  silicon  microelectronics[3,  4].  There  is  also 
interest  in  moving  the  lasing  from  the  far-  and  mid-infrared  range  to  the  near-infrared 
optical  communications  wavelengths,  A  =  1.3  or  1.55/i7n[5].  Since  the  latter  wavelength 
corresponds  to  a  photon  energy  of  SOOmeV,  the  Sii_xGex/Si  heterosystem  is  inadequate 
because  a  maximum  practical  valence  band  offset  of  only  ~  SOOmey  can  be  obtained 
for  pseudomorphic  Sii-xGe*  layers  at  x  =  0.5  ~  0.6.  Therefore,  alternative  large- 
bandgap,  nearly  lattice-matched  barrier  materials  for  Si  quantum  wells  must  be  sought; 
materials  with  sufficiently  large  band  offsets  with  respect  to  silicon.  Possible  candidates 
include,  ZnS,  BeSeTe,  CaFa,  SiO*,  SiOa,  the  Si/SiOa  superlattice,  and  7-AI2O3,  among 
others[5]-[7]. 

The  Si/ZnS  heterosystem  has  received  the  most  attention  as  current  advances 
in  epitaxy  technology  have  allowed  the  growth  of  heterostructures  consisting  of  polar 
and  nonpolar  materials  [8,  9].  The  lattice  mismatch  of  cubic  ZnS  with  respect  to  Si 
is  only  0.3%.  The  valence  band  offset  has  been  predicted  theoretically[10]-[13],  while 
recent  experiments[14]  show  that  the  value  is  close  to  1.5eV,  sufficiently  large  to  give 
intersubband  energy  differences  in  the  desired  SOOmey  range.  MBE  growth  of  ZnS 
upon  Si,  and  Si  upon  ZnS  have  been  demonstrated[9],  with  the  use  of  an  As  monolayer 
to  satisfy  the  local  bonding  requirements,  although  the  effect  of  the  monolayer  on  the 
offsets  has  not  been  determined. 

The  possibility  of  population  inversion  and  the  operation  of  intersubband  lasers 
depend  critically  on  the  lifetimes  of  the  involved  subbands.  The  subband  lifetimes  in 
turn  are  determined  by  nonradiative  phonon  scattering  processes.  The  purpose  of  the 
present  paper  is  to  study  the  optical  phonon  modes  and  their  interaction  with  carriers 
in  the  Si/ZnS  system  since  the  optical  phonon  scattering  is  considered  to  be  dominant 
in  the  phonon  scattering  processes.  This  combination  of  materials  is  new,  since  it  con¬ 
sists  of  both  a  nonpolar  and  polar  semiconductor.  Previous  studies  in  carrier  scattering 
by  confined  optical  phonons  in  heterostructures  have  been  focused  only  on  one  type  of 
phonons,  either  polar[15]-[23]  or  nonpolar [24]- [27].  In  the  current  situation  involving 
both  polar  and  nonpolar  materials,  carrier  scattering  by  both  types  of  phonons  needs 
to  be  considered.  To  the  best  of  our  knowledge,  there  has  not  been  any  reported  in¬ 
vestigation  on  this  mixed  nature  of  optical  phonons,  their  confinement  effect,  and  their 
interaction  with  carriers  in  a  heterostructure.  In  this  paper,  we  will  present  a  theoret¬ 
ical  study  based  on  the  macroscopic  continuum  model  to  describe  the  confined  optical 
phonon  modes  and  will  use  this  model  to  calculate  the  optical-phonon  scattering  of 
heavy  holes  in  a  heterostructure  consisting  of  polar  (ZnS)  and  nonpolar  materials  (Si), 
as  we  are  interested  in  the  feasibility  of  constructing  an  intersubband  laser  within  the 
valence  band  of  the  Si/ZnS  heterostructure.  This  valence  intersubband  laser  will  likely 
be  a  quantum-parallel  superlattice  laser[28]  or  a  quantum  cascade  laser.  Our  latest 
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thinking[29]  is  that  each  of  the  N  laser  periods  will  consists  of  one  square  Si  quantum 
well  containing  two  active  heavy-hole  subbands.  Due  to  the  non-parabolicity  of  these 
subbands,  we  believe  that  a  population  inversion,  localized  in  /r-space,  can  be  engineered 
between  the  subbands.  In  this  investigation,  we  will  consider  a  simple  superlattice  com¬ 
prised  of  alternating  layers  of  Si  and  ZnS,  much  like  the  flat-band  superlattice  of  the 
quantum  parallel  laser  (Fig.2(b)  of  [28]).  The  results  of  this  study  will  provide  the  basis 
for  the  calculation  of  subband  lifetimes  required  to  determine  laser  gain  and  threshold. 

As  described  below  in  greater  detail,  since  the  optical  dispersions  (frequency 
versus  wavevector)  of  the  silicon  {tuvsi  =  64meV')  and  zinc  sulphide  {hujzns  =  437716^) 
have  no  overlap,  the  optical  phonons  are  assumed  to  be  totally  confined  in  both  materials. 
In  the  silicon  layers,  a  continuum  model  with  double  hybridization  of  the  longitudinal 
optical  (LO)  and  transverse  optical  (TO)  modes  is  used  to  describe  the  vibration  patterns 
of  the  guided  modes[24].  The  only  boundary  condition  that  needs  to  be  satisfied  in  the 
Si  layers  is  the  vanishing  of  the  displacements  at  the  Si-ZnS  interface,  since  the  ZnS 
layers  can  be  considered  as  infinitely  rigid  with  respect  to  the  vibrations  of  the  Si  layer. 
Hence,  there  is  no  interface  mode  in  the  Si  layers.  The  situation  on  the  ZnS  layers  is  more 
complex.  Following  the  work  by  RidleyflG,  17],  here  a  continuum  model  is  employed 
with  hybridization  of  the  optical  LO,  TO,  and  interface  polariton  (IP)  modes  needed 
to  satisfy  both  the  mechanical  and  electrostatic  boundary  conditions  at  the  interfaces. 
Specifically,  the  electrostatic  boundary  conditions  are  the  continuity  of  Ei,  the  electric 
field  parallel  to  the  interface,  and  the  continuity  of  Dg,  the  displacement  field  normal  to 
the  interface.  The  mechanical  boundary  condition  is  again  the  vanishing  of  the  optical 
displacements  since  the  Si  layers  can  be  considered  as  infinitely  rigid  with  respect  to  the 
vibrations  of  the  ZnS  layers. 

Our  current  work  provides  a  complete  set  of  analytical  expressions  for  the  optical 
phonon  dispersion  relations,  optical  displacements,  and  associated  scalar  and  vector 
potentials.  These  expressions  are  subsequently  used  in  calculating  the  interaction  of 
heavy  holes  with  the  confined  optical  phonons.  In  Section  II,  we  establish  a  continuum 
model  for  the  optical  displacement  modes  in  Si  and  ZnS  layers  satisfying  both  mechanical 
and  electrical  boundary  conditions.  In  section  III,  we  outline  a  numerical  procedure  for 
determining  the  frequency  of  a  ZnS  optical  mode  inducing  the  intersubband  scattering. 
In  Section  IV,  we  describe  the  scalar  and  vector  potentials  associated  with  the  ZnS 
ionic  displacement  modes.  In  Section  V,  we  calculate  the  intersubband  scattering  rate 
due  to  the  emission  of  Si  and  ZnS  optical  phonons  since  the  emission  process  is  rather 
significant  compared  to  the  scattering  process  of  optical  phonon  absorption.  In  Section 
VI,  we  summarize  and  discuss  our  results  and  conclusions. 

II.  Mode  Patterns 

A  continuum  model  for  the  optical  modes  in  the  Si/ZnS  superlattice  is  employed. 
Both  mechanical  and  electrical  boundary  conditions  are  satisfied  at  the  heterointerfaces. 
Since  the  optical  dispersion  relations  (frequency  versus  phonon  wavevector)  in  the  two 
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bulk  materials  have  no  overlap,  the  phonons  are  taken  to  be  confined  in  their  respective 
materials.  For  the  Si  layers,  the  continuum  model  for  optical  phonons  in  nonpolar 
materials[24,  26]  is  used.  Here  double  hybridization  of  the  LO  (longitudinal  optical) 
and  TO  (transverse  optical)  modes  is  used  to  give  the  vibration  patterns  of  the  guided 
modes.  Since  the  ZnS  layers  are  infinitely  rigid  with  respect  to  the  vibrations  of  the  Si 
layers,  only  the  mechanical  boundary  condition,  the  vanishing  of  the  displacements  at 
the  interfaces,  has  to  be  satisfied. 

For  the  polar  ZnS  layers,  the  situation  is  more  complex  and  an  alternate  contin¬ 
uum  model[16,  17]  consisting  of  an  intermixing  of  confined  LO,  TO,  and  IP  (interface 
polariton)  modes  is  needed  in  order  to  satisfy  both  the  electrostatic  and  mechanical 
boundary  conditions.  The  boundary  conditions  which  must  be  satisfied  are  (1)  the  con¬ 
tinuity  of  Ex,  the  component  of  electric  field  parallel  to  the  interface,  (2)  the  continuity 
of  Dg,  the  component  of  the  displacement  vector  normal  to  the  interface,  and  (3)  the 
vanishing  of  the  vector  displacement  u  at  the  interface. 

A.  Modes  in  Si  Layers 

As  discussed  above,  since  the  ZnS  layers  can  be  treated  as  infinitely  rigid,  the 
boundary  condition  to  be  satisfied  in  the  Si  layers  is  the  vanishing  of  the  ionic  dis¬ 
placement  of  all  confined  vibration  modes.  This  is  an  assumption  of  strict  confinement 
yielding  only  the  guided  modes.  As  pointed  out  in  the  continuum  theory [24],  the  ionic 
displacement  of  confined  vibrations  has  two  components:  one  is  the  hybrid  of  the  LO 
and  p-polarized  TO  (p-TO)  modes,  and  other  is  the  uncoupled  s-polarized  TO  (s-TO) 
mode.  These  modes  are  defined  as  follows:  If  we  consider  a  {x,z)  plane  containing  the 
normal  to  the  layers  and  the  phonon  wavevector  Q,  then 

Q  ~  Qx^x  “b  Qz^z  (1) 

where  Cx  and  Cz  are  unit  vectors.  The  p-TO  mode  has  its  displacements  normal  to  Q 
and  in  the  plane,  while  the  s-TO  displacements  are  normal  to  Q  and  perpendicular  to 
the  plane  (||ey). 

The  form  of  the  ionic  displacement,  scalar,  and  vector  potentials  in  one  superlat¬ 
tice  period  differs  from  that  in  a  neighboring  period  only  by  a  phase  factor  proportional 
to  the  Bloch  superlattice  wavevector  qsl-  Their  expressions  given  below  are  obtained 
by  taking  qsl  =  0.  A  description  of  the  s-TO  mode  is 

Uy  =  +  B,.roe-‘«**),  (2) 

while  the  hybrid  of  the  LO  and  p-TO  modes  is  given  by 

Ux  =  [qx  ( Aioe’®''*  +  +  Qt  +  Bp-Toe~^’‘)] ,  ,o\ 

u,  =  -  Bioe"*’"")  -  gx(Vroe’’^"  -  Bp.Toe"'’^*)], 

which  axe  confined  within  the  Si  layer  with  a  width  oi  dsi,  0  <  z  <  dsi-  The  z- 
components  of  the  LO  and  TO  wavevector  have  been  distinguished  by  qi  and  qr,  re¬ 
spectively. 
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Since  the  LO  and  TO  modes  must  have  the  same  frequency  to  be  effectively 
coupled,  we  must  satisfy  the  condition 

u}^  =  Uq  —  =  u)q  —  oct{Qx  +  9t)>  (^) 


where  uq  is  the  bulk  Si  optical  phonon  frequency  at  F  point,  and  ar  are  the  sound 
velocities  of  LO  and  TO  dispersions  in  Si,  respectively. 

Using  the  boundary  condition  that  u  =  0  at  the  interfaces  gives  for  the  s-TO 

niT 

Uy  —  sm{qzz),  with  q^  —  —  (5) 

^Si 

where  n  =  1,2, ••  •  and  >1  is  a  mode  coefficient.  This  mode  does  not  mix  with  other 
modes. 


The  hybrid  LO  and  p-TO  modes  admit  two  classes  of  solutions. 

tion  is 

Ux  =  2Be^’^'^qx[cos{qiz)  -  cos(q'r2^)], 

o2 

Uz  =  sin(9i,2:)  +  —  sin  (972^)], 

Or 


The  ‘sine’  solu- 
(6) 


and  the  ‘cosine’  solution  is 


Ux  =  2iBe'^^^[qx  sm{qLz)  +  sin(97'2;)], 

9x 

Uz  =  2Be*®*®9i[cos(gi,2:)  -  cos(9t2)] 


where 


Ql  = 


nt'K 

dsi 


and 


TlxTT 


(8) 


where  ni,  =  1, 2,  •  •  ny  =  3, 4,  ■  •  tit  -  ni,  =  2, 4, 6,  •  •  *,  and  B  is  a  mode  coefficient. 
No  interface  modes  exist  in  the  Si  layer  because  of  the  boundary  condition  u  =  0. 


The  lowest  s-TO  mode  pattern  in  Eq.(5)  for  qz  =  'K/dsi  is  shown  in  Fig.  1(a) 
within  a  Si  layer  of  dsi  =  40A,  while  the  hybrid  patterns  of  the  lowest  p-TO  and  LO 
modes  with  q^  =  T^/dsi  and  qr  =  ^n/dsi  are  shown  in  Figs. 1(b)  and  1(c)  for  the  ‘sine’ 
and  ‘cosine’  solutions  given  in  Eqs.(6)  and  (7),  respectively  within  the  same  Si  layer.  The 
strict  confinement  which  requires  the  vanishing  of  ionic  displacements  at  the  boundaries 
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u  (arb.units) 


of  Si  layers  is  clearly  demonstrated  for  both  vibration  modes. 


Figure  1.  Vibration  patterns  in  a  Si  layer  with  a  width  of  40A  for  (a)  the  guided  s-TO 
mode,  (b)  the  ‘sine’  solution,  and  (c)  the  ‘cosine’  solution  of  the  guided  p-TO  and  LO 
modes. 


B.  Modes  in  ZnS  Layers 

The  boundary  conditions  are  the  continuity  of  Ex,  Dg,  and  the  vanishing  of  u  at 
the  interfaces.  These  conditions  can  be  satisfied  by  a  unique  linear  combination  of  LO, 
TO,  and  IP  modes  with  common  frequency  and  common  in-plane  wavevector  Qx, 


u  =  Uio  +  vlto  +  U/p.  (9) 

We  will  use  this  hybrid  expression  to  calculate  the  electrical  interaction  with  carriers 
which  is  considerably  stronger  than  the  optical  deformation  potential  interaction.  We 
need  consider  only  the  displacements  Ux  and  Ug,  since  Uy  associated  with  the  s-TO  mode 
has  no  related  electric  field  and  therefore  does  not  interact  with  carriers  electrically. 
Once  again,  the  expressions  are  obtained  by  taking  the  Bloch  superlattice  wavevector 
Qsl  —  0- 

For  the  LO  mode,  the  ionic  displacements  are 

Ux  =  -h  . 

Ug  =  ^  ^ 

which  is  confined  within  the  ZnS  layer  with  a  width  of  dzns,  — d^ns/2  <  z  <  dzus/"^ 
The  associated  electric  fields  are 


where 


Ex  —  Po'^X)  Ejg  —  Po^zy 


(11) 

(12) 
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with  the  effective  ionic  charge 


Coo 


where  M  is  the  reduced  mass,  €o  is  the  permittivity  of  free  space,  Coo,  e*  are  the  high- 
frequency  and  static  permittivities,  and  is  the  volume  of  primitive  unit  cell.  The  scalar 
potential  (f>  associated  with  the  electric  field  E  =  —  is  in  turn  given  as 


(14) 

For  the  TO  mode 

(15) 

The  electric  fields  associated  with  this  mode  are  negligible. 

For  the  IP  mode 

=  e‘(9*='-‘^*)gp(Ape‘9>’*  Bpe"*’*’^), 

(16) 

The  associated  electric  fields  are 

Bx  —  PpUi,  —  Ppf^zi 

(17) 

where 

_ 

Pp  —  Po  2  ,  ,2  > 

^LO  ~  ^TO 

(18) 

and  (Jlo  and  uro  are  bulk  ZnS  LO  and  TO  optical  phonon  frequencies  at  T  point, 
respectively.  The  electric  fields  associated  with  the  interface  modes  propagate  into  the 
Si  layers  although  they  are  treated  as  infinitely  rigid  and  do  not  contain  ZnS  ionic 
displacement. 

Being  a  transverse  electromagnetic  wave,  there  is  a  vector  potential  A  associated 
with  the  electric  field  E  =  —dA/dt.  Within  the  ZnS  layers, 

W)  (19) 

A,  = -  Bpe-‘«^*) 

UJ 

While  in  the  Si  layers,  a  similar  expression  can  be  obtained  with  another  set  of  mode 
coefficients,  Api  and  Bpi. 

Since  large  in-plane  wavevectors  are  likely  to  be  involved  when  dealing  with  carrier 
transitions  due  to  optical  phonons  between  two  subbands  separated  with  a  relatively 
large  energy,  we  have  endeavored  to  obtain  analytical  dispersions  of  the  LO  and  TO 
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optical  branches  by  curve-fitting  the  experimental  bulk  phonon  dispersions  for  the  entire 
Brillioun  zone.  The  requirement  for  common  frequency  yields, 

w  =  Wio  - /?l(9x  + 

=  uto  —  Pt{qI  + 

^  cjql  + 

n{uj) 

where  pi  =  0.808TFz  •  and  Pr  =  2.19^Hz  •  A^  are  obtained  from  curve-fitting 
the  bulk  ZnS  optical  phonon  dispersions,  c  is  the  velocity  of  light  in  vacuum,  and  n(a;) 
is  the  index  of  refraction.  In  the  above  expressions,  the  frequency  in  the  ZnS  layers 
lies  between  the  ZnS  LO  and  TO  zone  center  frequencies.  Since  loto  <  ^lo,  in  order 
for  the  TO  frequency  to  be  equal  to  a  LO  frequency  gr  must  be  imaginary  qr  =  iqg, 
corresponding  to  a  TO  interface  mode.  The  modes  which  interact  most  strongly  with 
carriers  are  those  with  frequencies  near  the  LO  branch.  For  these  modes,  the  value  of 
qo  is  large,  and  we  can  take  the  approximation 

tamhiqodzns)  «  1-  (21) 

In  the  unretarded  limit  (c  oo),  ql  +  qp^O  for  the  IP  mode.  Hence,  qp  w  iqx- 

Applying,  at  the  two  interfaces  between  layers  Si  and  ZnS  in  a  period  of  the 
superlattice,  the  conditions  that  Ux  and  Ug  equal  to  zero  along  with  the  continuity  of  Ex 
and  Dg,  leads  to  eight  simultaneous  equations  involving  the  eight  unknown  mode  coef¬ 
ficients  (Al,  Bl]  At,  Bt',  Ap,  Bp\  and  Api,Bpi).  The  following  two  ionic  displacement 
mode  patterns  emerge  for  the  hybrid  in  Eq.(9)  taking  the  Bloch  super  lattice  wavevector 
qsL  =  Q  and  the  approximation  taji'h{qodzns)  ^  1-  Both  ionic  displacement  patterns  are 
confined  within  the  ZnS  layer,  —dznsl^  <  z  <  dzns/^-  For  the  first  type. 


2iBe^^^^qx{sm{qiz) 

-[I -Pi  tanh(9xdzns/2)] siniqidzns/'^) 
2Be*^^^qL[cos{qLz) 

-  — Pi  sm{qLdzns/2)  — 

qi  '  cosh{qxdzns/2) 
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and  for  the  second  type, 

Ux  =  2Be*^^^qx{cos{qLz) 

-[l-p,  C0thferf2„5/2)] 

-P2  cos(5i.<J2„s/2)  5^77^^^}. 

smn{qxdznsl^) 

Uz  =  2zJ5e*®*®gj:,{sin(giz) 

+^[1  -  P2  coth(qxdzns/^)]  cos{qLdznsf‘2) smHqoz) 


,  Qx _ j  sinh(g^2)  , 


sinh(5odz„s/2) 


where 


Pi 

P2 

d 

r 


_  cosh{qxdsi/2)  cosh{qxdzns/2.) 

~  d 

_  smh{qxdsi/2)  smhjqxdzns/'^) 

d  ’ 

=  r  smh{qxdsi/2)  cosh{qxdzns/^)  +  sinh(gid^„s/2)  cosh(ga;dsi/2), 


and  epi  and  ep2  are  the  permittivities  in  Si  and  ZnS  layers,  respectively,  with 


Cp2  —  C(. 


-  ujIo 


(J-JPO 


(24) 


(25) 


To  illustrate  the  patterns  of  ionic  displacements  in  the  ZnS  layers  given  in  Eqs.(22) 
and  (23),  we  need  to  first  determine  values  for  qx,  qi,  and  qo-  To  do  so,  we  will  follow 
the  numerical  procedure  described  in  Section  III  by  arbitrarily  fixing  a  value  for  the  in¬ 
plane  phonon  wavevector  9*  =  'K/{hazns)j  where  azns  is  the  lattice  constant  of  ZnS.  In 
calculating  the  carrier-optical  phonon  interaction,  the  value  of  qx  is  actually  determined 
by  the  conservation  of  in-plane  momentum  between  the  initial  and  final  states  of  the 
scattering  process.  For  a  given  value  of  qx,  typically,  a  set  of  hybridized  modes  can  be 
obtained.  Here,  we  show  only  the  mode  pattern  with  frequency  close  to  ujio. 

We  obtained  huj  =  41meV',  qi  =  0.46  x  10®/cm  and  qo  =  0.48  x  10®/cm.  Substi¬ 
tuting  these  values  into  Eqs.(22)  and  (23),  we  obtained  Figs. 2(a)  and  2(b)  showing  the 
mode  patterns  of  ionic  displacement  of  both  the  first  and  second  types,  respectively  in 
a  ZnS  layer  of  dzns  =  20A.  It  can  be  seen  firom  Figs. 2(a)  and  2(b)  that  the  mechanical 
boundary  condition,  vanishing  of  the  ionic  displacements  at  the  interfaces  of  Si  and  ZnS 
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layers,  is  satisfied. 


z  (A)  Z  (A) 

Figure  2.  Vibration  patterns  in  a  ZnS  layer  with  a  width  of  20 A  for  (a)  the  first  type 
and  (b)  the  second  type  solutions  of  the  hybridized  LO,  TO  and  IP  modes. 

Ill.  Dispersion  Relationship 

The  phonon  frequency  in  the  ZnS  layers  is  determined  by  the  following  set  of 
equations; 

(  UJ  =  UJL0- MqI  +  qI), 

<  ui  =  ujto  -  Priql  -  Qo),  (26) 

^  ti+t^  C0s{qi,dzns)  +  ^3  S^^iqidzns)  =  0 


where 

ti  =  4psinh(gj;dsi)  +  4pr  sinhig^d  zns), 


h  =  —Apa, 

tz  =  Sp^r  smh{qxdzns)  siah{q^dsi)  -  4p^o? 

(27) 

+4p\‘^smh'^{qxdzns)  +  sinh^  {q^dsi)  +  1, 

with 

_ 

p  rrz  - 

Aqirsd 

(28) 

and 

^  )2  ,  ,2 
u/  ^TO 

^  9  9  * 

^LO  “  ^TO 

(29) 

The  third  equation  in  (26)  is  obtained  from  the  requirement  of  a  nonzero  solution  for 
the  eight  simultaneous  equations  discussed  above,  and  Eq.(27)  is  arrived  at  under  the 
approximation,  tanh(godzns)  «  1. 


The  numerical  procedure  for  determining  a  phonon  frequency  is  the  following: 
given  a  value  of  q^,  we  can  determine  those  of  ti,  <2,  and  from  Eq,(27).  Then  u  is 
scanned  from  ujto  to  ujio-  For  a  given  value  of  u,  qi  and  qo  are  obtained  from  the  first 


two  equations  in  (26).  Those  values  are  then  substituted  into  the  third  equation  in  (26) 
to  determine  if  the  particular  value  of  u;  is  a  solution. 

IV.  Scalar  and  Vector  Potentials 


Associated  with  the  two  types  of  ionic  displacement  in  Eqs.(22)  and  (23),  the 
scalar  potentials  in  ZnS  layers  are  given  as,  for  the  first  type, 


dzns 


ZnS  layer. 


2/9o5e*’*®sin(9L2i)  |2i|  <  ^ 

0  ^2!  <  %  Si  layer, 


(30) 


and  for  the  second  type. 


0  = 


dznS 


-2ipo5e’^*®cos(gi,2:i)  l-^il  <  — ZnS  layer, 


0 


22I  <  ^  Si  layer. 


(31) 


Note  that  we  have  used  two  different  coordinates  Zy  and  for  layers  ZnS  and  Si,  re¬ 
spectively,  with  their  origins  placed  at  the  centers  of  the  respective  layers. 


The  vector  potentials  can  be  obtained,  for  the  first  type, 

kl< 


u 

Mxpo 


'cosh(9idz„s/2) 


U) 

f  2ispoqx 
u 


sinh(ga;2;2) 


dznS 

k.i  <  f 


sin(gi,dzns/2) 


1. 

^^Be‘«-*Vicosh(9*Z2) 

u 


cosh(gxZi) 

cosh(9xdzns/2) 


1^1 1  < 


dznS 


ZnS  layer, 
Si  layer, 

ZnS  layer. 


(32) 


(33) 


and  for  the  second  type, 
2ispoqx 


Ax=< 


u 

^iQxPo 


cos{qLdzns/‘2) 


cosh{qxZi) 

smh{qxdzns/2) 


Zi\  < 


a; 

2spoqx 


cosh(9i2:2) 


Ax  = 


u 

MxPo 


J3e’«*®P2  cosiqidzns/^)- 


sinh(gx2i) 


sinh(9xd^„s/2) 


ki|< 


u 


Be*’*®V2  sinh(9iZ2) 


^  Si 

layer, 

dznS 

2 

ZnS  layer, 

dsi 

2 

Si  layer, 

( 

dznS 

2 

[ 

ZnS  layer, 

dsi 

2 

Si  layer. 

(34) 


(35) 
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where 


V,  = 


Vo  = 


rsin(gi,cizTv5/2)  cosh(gj:rfzns/2) 


r  COS 


(9rc?2„5/2^1iinh(gxrfzns/2) 


(36) 


2d 


and  pi,  P2,  and  d  are  given  in  Eq.(24). 


The  scalar  potentials  associated  with  the  LO  modes  are  strictly  confined  within 
the  ZnS  layers.  Their  distributions  are  shown  in  Fig.3  for  the  first  and  second  types 
given  in  Eqs.(30)  and  (31)  with  ql  =  0.31  x  10®/cm,  dznS  —  20A,  respectively. 

The  vector  potential  associated  with  the  IP  modes  are  distributed  in  both  Si  and 
ZnS  layers,  even  though  Si  layers  are  treated  as  infinitely  rigid  and  do  not  contain  ZnS 
ionic  displacements.  The  profiles  for  the  two  components  of  the  vector  potentials  given 
in  Eqs.  (32-35)  for  the  first  and  second  types  with  dsi  =  40A,  dznS  =  20A  are  shown  in 
Figs. 4(a)  and  4(b),  respectively. 


Figure  3.  Scalar  potential  distribution  associated  with  the  LO  modes  in  a  period  of 
the  Si/ZnS  superlattice  with  dsi  =  40 A  and  dznS  =  20 A  for  both  the  first  and  second 
types  of  the  vibration  modes. 

It  can  be  seen  from  Figs. 3  and  4  that  both  scalar  and  vector  potentials  are 
not  continuous  across  the  interfaces.  However,  as  pointed  by  Ridley  [16],  the  energy  of 
interaction  with  an  electron  traveling  coherently  with  the  optical  phonon  is  continuous. 
The  electric  field  can  be  obtained  as 

(37) 

The  continuity  of  Ex  and  Dg  =  e{u})Ez  implies  that  at  the  boundaries, 


^  =  aT- 


^^x\z2=±dsi/2  dZnS/2  |zi=^d2nS/2) 

^z\z2=±dsi/2  ~  '^^z\zi=:^dzn$ /2i 


(38) 
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where  Au  and  Au  axe  x-  and  z-components  of  the  vector  potential  in  Si  layers.  The 
interaction  in  the  Si  layer  is  e{AixVx  +  Auv^)  and  in  the  ZnS  layer  AxVx  +  AgVz), 

which  are  equal  when  the  electron  velocity  Vx  =  (jj/Qx  and  =  0.  Thus,  the  coherent 
interaction  energy  is  continuous  across  the  interfaces. 


Figure  4.  Vector  potentials  associated  with  the  IP  modes  distributed  in  a  period  of  the 
Si/ZnS  superlattice  with  dsi  =  40A  and  dzns  =  20A  for  (a)  the  first  type  and  (b)  the 
second  type  of  the  vibration  modes. 


The  electric  field  distributions  for  Ex  and  e(a;)  in  Si  {dsi  =  40A)  and  ZnS  {dznS  = 
40A)  layers  are  shown  in  Figs.5(a)  and  5(b)  for  the  first  and  second  types,  respectively. 
The  continuity  of  Ex  and  D*  across  the  Si  and  ZnS  interface  according  to  Eq.(38)  is 
clearly  demonstrated. 


Figure  5.  The  field  distributions.  Ex  and  Dg,  derived  from  the  scalar  and  vector 
potentials,  in  a  period  of  the  Si^ZnS  superlattice  ivith  dzi  40^1  and  dzns  ““  20A  for 
(a)  the  first  type  and  (b)  the  second  type  of  the  vibration  modes. 
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V.  Intersubband  Scattering 

Since  the  optical  modes  in  the  Si/ZnS  superlattice  consist  of  confined  nonpolar  Si 
and  polar  ZnS  optical  phonons,  the  calculation  of  carrier  scattering  by  optical  phonons 
in  such  a  structure  needs  to  include  contributions  from  both  types  of  phonons.  The 
Hamiltonian  that  describes  the  carrier  interaction  with  the  nonpolar  Si  optical  phonons 
is  given  by[30] 

ir=iD-u  (39) 

where  D  is  the  optical  deformation  potential.  This  Hamiltonian  obviously  vanishes 
outside  of  the  Si  layers  in  the  Si/ZnS  superlattice  since  the  Si  optical  displacement 
modes  are  strictly  confined  within  the  Si  layers.  However,  the  carrier  interaction  with 
the  confined  polar  ZnS  optical  phonons  extends  over  both  ZnS  and  Si  layers.  The 
electrical  interaction  Hamiltonian  can  be  obtained  using  the  scalar  and  vector  potentials 

H  =  -e<^  +  —  A  •  p,  (40) 

m 

where  p  is  the  momentum  operator,  e  and  m  are  the  free  electron  charge  and  mass, 
respectively.  Although  the  scalar  potential  (j)  associated  with  the  LO  mode  vanishes  in 
Si  layers,  the  A  •  p  interaction  exists  in  both  layers  since  the  vector  potential  associated 
with  the  interface  modes  in  the  ZnS  layers  as  shown  in  Fig.4  propagates  into  the  Si 
layers  as  well. 

A.  Scattering  due  to  Si  Phonons 


The  displacement  patterns  described  in  Eqs.(5-7)  all  contain  an  arbitrary  constant 
for  the  mode  amplitude  which  can  be  normalized  by  equating  the  energy  of  the  vibration 
mode  with  that  of  a  simple  harmonic  oscillator[16] 

where  S  is  the  sample  surface  area  (in  (x,  y)  plane),  Q  is  the  volume  of  the  unit  cell,  and 
X  is  the  normal  coordinator  of  the  oscillator.  The  heavy-hole  state  can  be  characterized 
by  |k,  n  >  with  the  in-plane  momentum  k  and  subband  index  n.  In  the  approximation 
of  constant  eflfective  mass  for  heavy  holes,  the  matrix  element  for  the  transition  from 
state  |k,  n  >  to  |k',  n'  >  due  to  the  emission  of  a  nonpolar  Si  optical  phonon  is 


<  k',n'|iflk,n>=  { 


2ps!uoSdliAA\q,)  (?*)  (s-TO) 

I  ninjup)  +  1]  ^ 

‘^Psi<^oSdsiAc{qLi  Qt) 

■  l^xGln'  {QL,9r)  +  Dz  (qi ,  qr)]  (hybrid) , 


(42) 


for  the  s-TO  mode  and  the  hybrid  of  the  LO  and  p-TO  mode,  respectively.  n{uo) 
is  the  number  of  Si  optical  phonons  at  thermal  equilibrium,  and  psi  is  the  density 
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of  Si.  The  three  components  of  the  optical  deformation  potential,  Dy,  and  Dg 
are  assumed  equal  to  =  D/\/3  in  the  calculation,  in  view  of  the  assumption  of 
isotropy.  The  Kronecker  symbol  indicates  the  in-plane  {x,  y)  momentum  conservation. 
The  normalization  factors  are  given  by 


^A{qz)  =  -^  f  ^'u*yUydz  (s-TO), 

Ac (9l I  +  u].u^)dz  (hybrid) . 

asi  JQ 


(43) 


The  Gnn'-fuiictions  contain  envelope  wavefunctions,  ipn  and  from  which  interference 
effect  can  be  obtained.  Specifically, 


(44) 


for  the  s-TO  mode,  and 


GL' (9n.fr)  =  [  ^'^n^n'Uxdz, 

(fr.fr)  =  /  ^''iPnfpn'Uzdz, 

J  0 


(45) 


for  the  hybrid  of  LO  and  p-TO  mode.  The  heavy-hole  energy  levels  and  envelope  wave- 
functions  are  obtained  by  the  finite  square  well  model  for  the  superlattice  with  the 
heavy-hole  band  offeet  taken  to  be  1.5eV'[14]. 


Applying  the  Fermi  golden  rule,  we  obtain  the  scattering  rate  due  to  the  emission 
of  a  nonpolar  Si  optical  phonon. 


f  ^]-^o  Y'  Idlin' I  (c  TOl 

.  2n^p,Uodsi  A^  ^  ’ 


(46) 


where  we  have  assumed  that  for  the  intersubband  process  (n  ^  n')  the  heavy  holes  are 
scattered  from  the  bottdm  of  their  original  subbands  (k  =  0),  and  the  sum  is  over  those 
participating  modes  of  Eq.(8)  that,  according  to  Eq.(4),  yield  values  of  qx  satisfying  the 
in-plane  momentum  conservation [24]. 


B.  Scattering  due  to  ZnS  Phonons 


The  normalization  of  the  amplitudes  of  the  confined  ZnS  displacement  modes 
can  be  carried  out  by  equating  the  energy  of  a  hybrid,  a  mixture  of  mechanical  and 
eletromagnetic  energies,  with  that  of  a  simple  harmonic  oscillator[16].  Since  only  the  IP 
mode  contributes  electromagnetic  energy  which  is  small  in  magnitude  when  compared 
with  the  mechanical  energy,  neglecting  it  entirely  will  introduce  little  error  in  evaluating 
the  energy  of  a  hybrid  ZnS  mode.  We  therefore  can  use  Eq.(41)  for  the  normalization 
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of  a  ZnS  mode,  except  that  now  the  integral  is  over  the  ZnS  layer,  dzns-  The  matrix 
element  for  the  transition  from  state  |k,  n  >  to  |k',n'  >  due  to  the  emission  of  a  polar 
ZnS  optical  phonon  is 


<k',n'lJ?|k,n>= 


\ 


e  h.[7l(w2ns)  “t"  l-]^Zn5  j-  (  z^tx  /^z  \  /a'7\ 


for  both  the  first  and  second  types.  n{uizns)  is  the  number  of  ZnS  optical  phonons  at 
thermal  equilibrium,  and 

i  =  — -i.  (48) 

Cp  Coo 

The  normalization  factors  for  both  the  first  and  second  types  can  all  be  calculated  by 


1  fdznS 

Ai,2  =  -j —  /  [ulu^  +  u*u*)d2r 
dZnS  Jo 


(49) 


with  optical  displacements  given  in  Eqs.(22)  and  (23).  The  Gnn' -functions  containing  the 
interference  effect  between  two  subband  envelope  wavefunction,  and  V’n')  given 
specifically  as 

<^nn'  =  i-  /  4>Tp^1pn’dz,  (50) 

Pa 


for  the  scalar  potential  scattering  associated  with  the  LO  modes,  and 

r^ZnS 
/O, 


1  raZnS 

Gin'  =  A^iPliPn'dz, 

1  fdznS  d 

Gin'  ~  J^z'>Pl-^‘fpn'dz, 


(51) 


for  the  vector  potential  scattering  associated  with  the  IP  modes.  Applying  the  Fermi 
golden  rule,  the  intersubband  scattering  rate  due  to  the  emission  of  a  polar  ZnS  phonon 
can  then  be  obtained  by  taking  summation  over  contributing  confined  ZnS  optical  modes 


_  e’^mlf^lnjujLo)  -F  Ijo^LO 
^^^dznS^p 


G0  12 

Ti.Ti'  I  ^nn'  _ I 


m 


m 


Ai,2 


(52) 


where  we  have  again  assumed  that  the  heavy  holes  are  scattered  from  the  bottom  of 
subband  n  (k  =  0),  and  have  taken  the  approximation  of  uzns  =  ojio  since  the  modes 
which  interact  most  strongly  with  carriers  are  those  with  frequencies  near  the  LO  branch. 


C.  Intersubband  Scattering  Rates 


The  scattering  rates  due  to  the  emission  of  Si  and  ZnS  optical  phonons  were 
calculated  for  the  intersuband  transition  (2-1)  originated  from  the  bottom  of  the  heavy- 
hole  subband  2  with  zero  kinetic  energy  to  heavy-hole  subband  1.  Figure  6  shows  the 
2-1  scattering  rates  as  a  function  of  the  Si  well  width  while  fixing  the  barrier  width  at 
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dzns  =  40A  in  the  Si/ZnS  superlattice.  The  total  scattering  rate  is  the  summation  of 
contributions  from  the  heavy-hole  interaction  with  Si  and  ZnS  optical  phonons.  In  the 
small  well  width  region  {dsi  <  30A),  the  heavy-hole  scattering  due  to  the  ZnS  optical 
phonons  is  stronger  than  that  due  to  the  Si  optical  phonons.  This  is  attributed  to  the  faet 
that  when  the  Si  well  width  is  small  there  is  significant  envelope  function  overlap  between 
subbands  1  and  2  in  the  ZnS  barrier  region  where  the  ZnS  LO  phonons  are  confined.  As 
the  well  width  increases,  the  distribution  of  envelope  functions  in  the  barrier  decreases. 
As  a  result,  the  scattering  due  to  the  ZnS  LO  phonons  reduces  considerably  and  the  ZnS 
phonon  scattering  is  mostly  through  IP  modes  which  propagate  throughout  the  super¬ 
lattice  structure.  As  the  well  width  continues  to  increase,  the  energy  separation  between 
subband  1  and  2  decreases.  The  intersubband  scattering  between  these  two  subbands 
requires  an  emitted  ZnS  phonon  with  a  small  in-plane  wavevector  (qx)  in  order  to  satisfy 
the  in-plane  momentum  conservation  for  the  scattering  process  to  take  place.  This  leads 
to  an  increased  intersubband  scattering  rate  since  polar  optical  phonons  with  smaller 
wavevectors  interact  more  strongly  with  carriers  to  induce  intersubband  transitions  as 
suggested  by  the  well-known  l/(g^  -f  dependence  of  the  interaction  Hamiltonian  in 
polar  material  quantum  wells[15].  A  similar  dependence  of  the  intersubband  scattering 
rate  of  Eq.(52)  due  to  the  confined  ZnS  optical  phonons  is  implicitly  included  in  the  nor¬ 
malization  factor  (A  1,2)  given  by  Eq.(49).  Further  increasing  the  well  width  to  dsi  >  82 A 
causes  the  energy  separation  between  the  two  subbands  to  be  less  than  the  ZnS  optical 
phonon  energy  (4377161^)  and  the  heavy  holes  at  the  bottom  of  subband  2  cannot  emit 
ZnS  optical  phonons  to  make  a  transition  to  subband  1 ,  resulting  in  zero  scattering  rate 
due  to  the  emission  of  ZnS  optical  phonon.  The  scattering  rate  due  to  the  emission 
of  Si  optical  phonon  confined  within  the  Si  well  demonstrates  a  steady  decrease  as  the 
well  width  (d^j)  increases,  which  suggests  that  the  factor  1/dsi  in  Eq.(46)  dominates 
the  small  increase  in  the  interference  effect  G„n'  function.  As  the  well  width  increases 
beyond  62A,  the  energy  separation  between  the  two  heavy-hole  subbands  becomes  less 
than  the  Si  optical  phonon  energy  (64TneV').  As  a  result,  the  scattering  rate  due  to  the 
emission  of  Si  optical  phonons  reduces  to  zero,  in  which  case  the  heavy-hole  lifetime 
of  subband  2  can  be  enhanced  dramatically  since  the  significant  scattering  process  of 
optical  phonon  emission  is  suppressed  although  the  weaker  optical  phonon  absorption 
and  acoustic  phonon  scattering  processes  are  still  possible. 

Figure  7  shows  the  intersubband  scattering  rates  between  the  same  two  heavy- 
hole  subbands  due  to  the  emission  of  both  Si  and  ZnS  optical  phonons  as  a  function  of 
the  barrier  width  (dzns)  in  the  Si/ZnS  supperlattice.  The  well  width  (dsi)  is  fixed  at 
30A.  The  scattering  rate  due  to  the  emission  of  Si  optical  phonons  remains  unchanged 
as  the  barrier  width  varies  since  the  subband  energy  levels  are  hardly  shifted  and  the 
G„n»-function  for  the  Si  phonon  scattering  has  little  noticeable  change.  The  scattering 
rate  due  to  the  emission  of  an  ZnS  optical  phonon,  on  the  other  hand,  demonstrates 
a  decreasing  trend  as  the  ZnS  barrier  width  increases  as  suggested  in  Eq.(52)  with  the 
factor  of  1/dznS-  The  small  discontinuous  incremental  steps  in  the  ZnS-scattering  curve 
are  due  to  the  discrete  nature  of  the  increase  in  the  number  of  allowed  LO  modes  confined 
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in  the  ZnS  barrier  as  it  increases. 


Si  layer  thickness  (A) 

Figure  6.  Intersubband  scattering  rates  due  to  the  emission  of  Si  and  ZnS  optical 
phonons  as  a  function  of  Si  well  width  (dsi)  for  a  barrier  width  of  dzns  =  40 A. 


Figure  7.  Intersubband  scattering  rates  due  to  the  emission  of  Si  and  ZnS  optical 
phonons  as  a  function  of  ZnS  barrier  width  idzns)  for  a  well  width  of  dsi  —  30A. 

VI.  Summary  and  Discussion 

We  have  provided  an  analytical  model  of  optical  modes  in  Si/ZnS  superlattices 
consisting  of  polar  and  nonpolar  optical  phonons.  This  is  a  new  procedure  for  obtaining 
the  eigen  modes  of  a  mixed  polar-nonpolar  heterosystem.  In  the  Si  layers,  a  continuum 
model  with  double  hybridization  of  the  LO  and  TO  modes  is  used  to  describe  the 
vibration  patterns.  Since  there  is  no  electric  field  resulting  from  the  nonpolar  ionic 
displacements  in  Si  layers,  the  only  boundary  condition  that  needs  to  be  satisfied  in  the 
Si  layers  is  the  vanishing  of  the  displacements  at  the  Si-ZnS  interface,  as  the  ZnS  layers 
can  be  considered  as  infinitely  rigid  with  respect  to  the  vibrations  of  the  Si  layer.  Due  to 
this  strict  confinement,  only  guided  modes  emerge  in  the  Si  layers  which  consist  of  s-TO 
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and  coupled  p-TO  and  LO  modes,  with  no  interface  modes.  These  guided  modes  have 
been  illustrated.  Their  interaction  with  carriers  in  the  superlattice  can  be  calculated 
through  the  optical  deformation  potential  for  Si.  The  interaction  Hamiltonian  can  be 
obtained  by  talcing  the  product  of  this  potential  with  the  normalized  ionic  displacement. 

However,  for  the  optical  phonons  in  ZnS  layers,  we  need  to  include  the  electrical 
interaction  in  calculating  the  carrier  scattering  by  optical  phonons,  since  there  are  elec¬ 
tric  fields  associated  with  the  polar  optical  vibrations.  As  a  result,  both  mechanical  and 
electrostatic  boundary  conditions  need  to  be  satisfied  in  the  interfaces.  A  continuum 
model  employing  a  linear  combination  of  LO,  TO  and  IP  (interface  polariton)  modes 
with  a  common  frequency  is  used  to  describe  the  ionic  displacements  in  ZnS  layers.  A 
numerical  procedure  for  determining  a  phonon  frequency  is  provided.  This  hybridized 
model  is  necessary  to  meet  the  simultaneous  requirement  on  the  mechanical  and  electro¬ 
static  boundary  conditions.  The  mechanical  boundary  condition  is  again  the  vanishing 
of  the  optical  displacements  since  Si  layers  can  be  considered  as  infinitely  rigid  with 
respect  to  the  vibrations  of  the  ZnS  layers.  The  electrostatic  boundary  conditions  are 
the  continuity  of  the  electric  field  parallel  to  the  interface,  and  the  continuity  of  the 
displacement  field  normal  to  the  interface.  Based  on  this  set  of  boundary  conditions, 
expressions  are  obtained  for  the  ionic  displacements  in  ZnS  layers  consisting  of  LO,  TO, 
and  IP  modes.  There  are  scalar  and  vector  potentials  associated  with  the  LO  and  IP 
modes,  respectively,  but  no  electric  field  associated  with  the  TO  mode.  The  scalar  po¬ 
tential  and  its  associated  electric  field  due  to  the  LO  mode  are  distributed  only  within 
the  ZnS  layers  and  are  zero  in  the  Si  layers.  But  the  vector  potential  and  its  associ¬ 
ated  electric  field  due  to  the  IP  mode  have  distributions  in  both  ZnS  and  Si  layers  even 
though  there  is  no  ZnS  ionic  displacement  mode  in  the  Si  layers.  Examples  of  these 
mode  characteristics  have  been  demonstrated.  Neither  the  scalar  nor  vector  potential 
is  continuous  across  the  Si-ZnS  interface,  but  the  energy  of  coherent  interaction  with 
carriers  is  continuous  due  to  the  continuity  of  the  electric  field  parallel  to  the  interface. 

The  analytical  model  for  the  confined  optical  modes  consisting  of  polar  and  non¬ 
polar  optical  phonons  is  employed  in  calculating  the  carrier-phonon  interaction.  Our 
results  indicate  that  contributions  to  heavy-hole  intersubband  scattering  from  confined 
Si  and  ZnS  optical  phonons  strongly  depend  on  the  well  width  since  it  varies  the  dis¬ 
tributions  of  envelope  functions  of  involved  subbands  which  ultimately  determines  the 
intersubband  scattering  between  them  through  the  overlapping  interference  effect  Gnn- 
function.  For  small  Si  well  width  (<  30 A),  the  scattering  rate  due  to  ZnS  optical  phonon 
is  stronger  than  that  of  Si  optical  phonons.  As  the  well  width  increases  the  scattering 
rate  due  to  the  Si  optical  phonons  surpasses  that  of  ZnS  optical  phonons.  The  scattering 
rate  dependence  on  barrier  width  is  relatively  weak. 
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